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Abstract 

The degenerate Fermi gas coupled to a random potential is used to 
study metal-insulator transitions in various dimensions. We first recast 
the problem in the sea-boson language that allows for an easy evaluation 
of important physical attributes. We evaluate the dynamical number- 
number correlation function and from this compute the a.c. conductivity. 
We find that the d.c. conductivity vanishes in one and two dimensions. 
For a hamiltonian that forbids scattering of an electron from within the 
Fermi surface to another state within the Fermi surface we find that there 
is no metal-insulator transition in three dimensions either. 



1 Introduction 

In a series of published works [1] [2], and a recent preprint [3], we showed how 
to extract the anomalous exponents in case of the Luttinger model using sea- 
bosons. This paves the way for application of the amended sea-boson theory 
that now is powerful enough to reproduce most of the exactly known results in 
Id, to other systems such as electrons with quenched disorder with and with- 
out Coulomb interactions in various dimensions. The relevant literature on this 
subject is vast and we shall not attempt to be exhaustive in surveying it. Ander- 
son's pioneering work on localization [4] was followed by the work of Abrahams 
et.al. [8] and later on a more rigorous formulation of the notion of disorder av- 
eraging was given by McKane and Stone [6]. This relates to a single electron in 
a disordered potential. The classic review of Lee and Ramakrishnan [7] includes 
many references on the literature concerning the degenerate electron gas in a 
disordered potential. A more recent review is by Abrahams et. al. [8]. 
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2 Number-Number Correlation Function 



Eventually, we would like to compute the a.c. conductivity at absolute zero. 
Unfortunately this quantity is rather difficult to compute. This is because it 
involves first calculating the dynamical number-number correlation function. 
This latter function has proved very difficult to evaluate. Before we evaluate 
this quantity we would like to say a few words about how the dynamical number- 
number correlation enters into the picture. It is defined as follows. 

N(kt;k't) =< n k (t)n k ,(t) > - < n k (t) >< n k ,(t) > (1) 

Notice that the a.c. conductivity is related to the dynamical total momentum- 
momentum correlation function. This formula was derived in an earlier preprint [3] . 
The momentum-momentum correlation function in turn may be related to the 
dynamical number-number correlation function. 

(<SP(t)- 5P(0))=£(k.k') [(n k (f)n k ,(0)}-(%(t))K'(0)>] (2) 
kk' 

In a recent preprint [3], wc provided some hints as to how might go about com- 
puting the number-number correlation function for the interacting system. It 
involves functional differentiation of the average momentum distribution with 
respect to sources that couple to the number operator. When this is done care- 
fully we find the following appealing form of the dynamical number-number 
correlation function. 

N(kt; kY) = (1 -n F (k))(l - n F (k'))S A AQst; kY) + n F (k)n F (k )S B B(kt; kY) 

-(1 - n F (k))n F (k)S A B(kt; kY) - (1 - n F {k ))n F (k)S B A(kt; kY) (3) 
Here the various quantities are defined recursively. First (here m,n = A, B), 

S mn (kt;k't) = e -2<s m (k)> e -2<5„(k')> S ° mn (kt;kt) (4) 

One could take the point of view that S , r °„„(kt;kt ) is evaluated by assum- 
ing that otk(q) are canonical bosons, dropping all the square roots and so on. 
The reason being that the corrections caused by fluctuations in the momentum 
distributions are included in the exponential prefactors. These quantities are 
defined recursively. 

S A (k, t) = J2 4-q/2 (<*> *)^k-q/2 (q, *) (5) 

q 

S B (k, t)=J2 4+q/2 (q ' f ) A k+q/2 (q, t) (6) 

q 

S° mn (kt;kt) = (s m (k,t)S„(k',t')) - (S m (k,t)) (S n (k,t)) (7) 
where m,n — A, B. 
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3 The Toy Hamiltonian 



Here we couple the free Fermi gas to a disorder potential and compute the a.c. 
conductivity. The diagonalization is rendered trivial in the sea-boson language. 
However, the formula for the dynamical number-number correlation function in 
terms of the bosons is very nontrivial and can therefore be expected to lead to 
nontrivial results. 



H = H ^4(q)^(q) + E ^# E + 4(q)] (8) 

kq n V i. 



The above hamiltonian describes electrons close to the Fermi surface interacting 
with the disorder potential. However, notice that no externally chosen cutoff is 
needed. A natural smooth cutoff emerges by not linearizing the bare fermion 
dispersion. In the Fermi language, Eq.( 8) is equivalent to the following hamil- 
tonian. 

H = E ek4c k + £ [A k (-q) + A k (q)] Ck+q/2 c k - q/2 (9) 

k k,q 

Here A k (q) = n F (k + q/2)(l - n F (k- q/2)) and n F (k) = 0(fc F -|k|). Thus the 
toy hamiltonian Eq.( 9) describes electrons coupling to the disorder potential 
near the Fermi surface in such a way that processes that take an electron below 
the Fermi surface and place it in another state also below the Fermi surface 
or both above the Fermi surface are forbidden. We shall see that in this case 
there is no metal insulator transition in any dimension. However, we reproduce 
the results that in one and two dimensions, the d.c. conductivity is zero. This 
hamiltonian may be trivially diagonalized by the following transformation. 

Ak(q,*)=A k (q)e-^-^^ (10) 



Thus we have, 

V (k.q 

Also for the number fluctuations 



4(q)^(q)> = ' T ;| nrrW * - q/2)(i - Mk + q/2)) (11) 



4(q, t)A k (q, t)A[, (q, t')A k ' (q, *'))-( A^q, t)A k (q, t)) U], (q , t)A^ (q , t) 



= '^ q) ' (k m q)2 e-^(*-*')5 kik ^ qiq m F (k-q/2)(l-n F (k + q/2)) (12) 
Thus we may compute the following quantities, 

< &(k) >= E |t/rf y q)|2 (( k -q/ 2 ),q) 2 nf(k - q)(1 - nF{k)) (13) 
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< §bM >= E lUd i q) \ (k+ ^ /2 y qr -^ - + q)) (14) 
l^(q)| 2 1 i( >La_ eq )( t _ t ' ) 



S AA (k,t;k',i') = 5 k;k , E T , kq _ , 2 e-^-^ > nF (k-q)(l-n,(k)) 
q v I— £ qi 

(15) 

S BS (k,t;k\i') = <5 k;k , E '^ S i q)l \ k . q 1 , 2 e-^^)(*"*') „,(k)(l-n,(k+q)) 

q I m + £ q) 

(16) 

SA B (k,t;k',t') = l £ W k - k ')l 2 1 e -^-v)(*-*') nF (k')(l-n F (k)) 

(17) 

S Bjl (M;k',t') = l^"( k ~ k )l 2 ^— e^V-^)(*-*')^(k)(l-n F (k')) 



(18) 

In an earlier preprint we showed that the real part of the a.c. conductivity may 
be written as, 

Re U dls )] = (^f) ^D tk ') N(k,e tl e r ,k ,0) S(uj-e t + ej ) (19) 

kk' hi 

where, 

JV(k,i;k',0) = (n k (t)n k /(0))-<n k (t)) (n k -(0)) = ^ ^^^(k, e,, ; k' , 0) 

(20) 



g AA (k,t;k,t , ) = ^ kk ,E^, k |t/rfis(k I / kj)l \ , 2 e-^-^)(*- t ')n F (k 3 )(l-n f (k i )) 

(21) 



&u,(M;k',t ) = E^,k'^, k '^ (k ; kj) ' 2 , . , 2 e-^^)(^)n f (k,)(l-n f (k0) 



(22) 

V (e.-e,-) 2 



Wk,t;k , ,f , ) = X)^,k'^^ l ^" (kt y kj) ' 2 r ,2 e "' (£ '" £j )(t - t ' ) ^(k J -)(i-n F (k I )) 



(23) 



S BB (k,i;k',*') = * k ,k'E^ V (e-V/ ^"" W(k,)(l-n,(kO) 

(24) 



jV(k, £t , ej ;k',0)^ e -^( k )^ k ^ ki , k l ^ (k ; kj)|2 1 n f (k 3 )(l-n f (k,)) 

V (Cj Cj) 

_ e -25i( k ) e - 2 50 (k ')^ k/ ^ ik l^(k.-k 3 .)| 2 1 F(k , )(1 _ nF(kj)) 

-r^« k »r-V^ ki , 4j , k M_^ 1 (kj)(1 _ „ F(ki)) 
+e-^( k ) (5k;k< ^. ;k i^( k ^-^)l 2 1 (k . )(1 _ nF(ki)) (25) 

4 A.C. Conductivity 

The disorder averaged a.c. conductivity for Gaussian disorder may be writtten 

as 

=< Re [<j(uj; U dis )} > dls 

= (^) ^E( k ? - isi(t " +*S -«-' k ") ^tt^mmu-m*.)) ^-S+S> 

(26) 

where, 

= E ^ (k . q m q2 /2)2 ^(k - q)(l - *F(k)) (27) 

^( k ) = £ ^ ( k . q rqV2) 2 nF(k)(1 - ^ (k + q)) (28) 
In other words, 

<j(uj) =< Re [<j(w; U dis )} > dls 

= (^f)h[ F **><r+w<) ((-+ e ) e - 4 ^ +£) +-- 4 ^ (£) ) 

(29) 

S^(e) = A 2 f F deD{e) (30) 
S°(e) = A 2 ^ deD(e)j-r—^ (31) 



^ (e)= (^F (2me) ^ (32) 

Using Mathematica™ we find, 

S°( £ ) = A 2 ^(2m)^ (^pj H2Fl[2,d/2,l + d/2, e -f] (33) 

( _n (-4+d)/2\ 

S°( e ) = A 2 ^(2m)^ f ( ^ + d) j H2Fl[2,2-d/2,3-d/2,^-] (34) 
In one dimension, further simplification is not possible. Hence we write, 

<$( c ) = A 2 ^(2m)-*^^ H2F1[2, 1/2,3/2,^] (35) 

S°(e) = A 2 ^(2m)-*^f^ H2Fl[2,3/2,5/2,^-] (36) 
In two spatial dimensions we have, 

"^(^"l) ,37) 

In three spatial dimensions we have, 

^( £ ) = A 2 ^3(2m)^ (^pj ff2Fl[2, 3/2, 5/2,^] (39) 

S£(e) = A 2 ^(2m)* ( 2£f 1/2 ) iT2Fl[2, 1/2, 3/2, ±] (40) 
In two dimensions we have, 



2^e 2 A 2 \ J_jn 2 f' F 
m J uj 3 (2 7 , , 



(41) 



This may be approximately evaluated as follows. 

, . /27re 2 A 2 \ 1 to 2 / _ 4 a 2 -££ t £\ 

It can bee seen that the zero frequency limit of the above expression is zero since 
the integral vanishes exponentially fast ~ e~ c °l w jij 1 . Thus the d.c. conductivity 
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of a two dimensional system is zero and the frequency dependence is rather 
nontrivial. Similarly we may expect that in one dimension the d.c. conductivity 
vanishes. Unfortunately for a similar reason we find that the d.c. conductivity 
in three dimensions also vanishes. This means we have to include terms beyond 
what Eq.( 9) does. Perhaps the reader can do this or at least offer to collaborate 
with the author. Please contact me at gsetlur@imsc.res.in 

5 Some Technical Musings 

It appears that the mathematical literature on the subject of quantum par- 
ticles in random potentials is vast [9]. It is possible, indeed likely that many 
mathematically rigorous results are known regarding this problem. But this 
does not prevent the authors from making some remarks that more knowledge- 
able readers may choose to critique. In particular, the author is uncomfortable 
with the notion of disorder averaging. Nature chooses its potentials based on 
the distribution of impurities, defects and so on. This potential is fixed and 
well-defined for a particular distribution of these imperfections. The physicists' 
ignorance of the precise nature of this potential is not a license to average over 
these potentials. Nature does not average, people do. But are people justified in 
averaging ? In other words can averaging simplify the problem without washing 
out essential physics ? In order to answer this question we have to make the 
following conjectures. 

DefnO : Let Ud be the set of all potentials U(x) in a fixed spatial dimension d. 
Defnl : Let Td be the set of all potentials C/(x) in a fixed spatial dimension d 
that has the following property. They all lead to the same exponent 5 for the 
frequency dependence of the a.c. conductivity. In other words, each of these 
potentials predicits that Re[a(u>)] ~ uj s (in some region of u> with possibly some 
additive part independent of u) with the same 5. 
Conjecture 1 : Td is dense in Ud- 

If Conjecturel is valid, then one may average over all these 'sufficiently er- 
ratic' potentials and expect to extract 5 which is all that physicists care about. It 
is possible that S may be extracted from a numerical solution of the Schrodinger 
equation using a specific U that belongs to the set Td- But this would involve 
using the computer for more than checking one's email, and not everyone likes 
that. 

Defn2 : Let Ms be the set of all potentials U(x) in spatial dimension d = 3 
that has the following property They all lead to the same exponent /3 for the 
mobility edge exponent. In other words, each of these potentials predicit that 
&d.c. ~ {Ef — E c )P 6{E f — E c ) with the same [3. However for different potentials, 
E c - the mobility edge, may be different. 
Conjecture 2 : JF 3 is dense in U3. 

If Conjecture2 is valid, then one may average over all these 'sufficiently 
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erratic' potentials and expect to extract /3. 

Thus the validity of the process of averaging over potentials rests crucially 
it seems, on all these sufficiently erratic potentials predicting the same expo- 
nents and on these sufficiently erratic potentials spanning nearly all possible 
potentials. 

If both these are satistifed then one may average over all potentials and 
extract the exponents, or, if one is better at programming, choose a particular 
potential from this set, numerically solve the Schrodinger equation and extract 
the exponent from there. In either case we should get the same answer. A final 
conjecture seems appropriate. 

Conjecture 3 : Let M 3 have an exponent /3 and T d have an exponent S , 
then (3 — (3 and 5 — 5 . In other words, these exponents are unique. 

With powerful computers now available, purely analytical methods such as 
this work may seem passe, but a closed formula for the a.c. conductivity that 
one can stare at (and one that is hopefully right) and admire has a charm that 
a cold data file on the hard disk is unable to duplicate. Besides, with Coulomb 
interaction, the problem becomes intractable numerically, however, one may 
expect to combine the sea-boson method with the present one to extract the 
exponents analytically. 

References 

[1] G.S. Sctlur and Y.C. Chang, Phys. Rev. B, 57, 15144 (1998). 
[2] G.S. Setlur and D.S. Citrin, Phys. Rev. B, 65, 165111 (2002). 
[3] G.S. Sctlur, cond-mat/0209610, cond-mat/0305039. 
[4] P. W. Anderson, Phys. Rev., 109, 1492 (1958). 

[5] E. Abrahams, P.W. Anderson, D. C. Licciardcllo and T.V. Ramakrishnan, 
Phys. Rev. Lett. 42, 673 (1979). 

[6] A. J. McKane and M. Stone, Ann. Phys. (N.Y.) 131, 36 (1981). 

[7] P.A. Lee and T.V. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985). 

[8] E. Abrahams, S.V. Kravchenko, M.P. Sarachik, Rev. Mod. Phys., 73, 251, 
(2001). 

[9] M. Krishna and K. B. Sinha, Proc. Ind. Acad. Sci. (Math. Sci.) Ill, 179 
(2001). 



8 



